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We put forward a new model-independent reconstruction scheme for dark energy which utilises the 
expected geometrical features of the luminosity-distance relation. The important advantage of this 
scheme is that it does not assume explicit ansatzes for cosmological parameters but only some very 
general cosmological properties via the geometrical features of the reconstructed luminosity-distance 
relation. Using the recently released supernovae data by the Supernova Legacy Survey together with 
a phase space representation, we show that the reconstructed luminosity-distance curves best fitting 
the data correspond to a slightly varying dark energy density with the Universe expanding slightly 
slower than the KCDM model. However, the ACDM model fits the data at la significance level 
and the fact that our best fitting luminosity-distance curve is lower than that of the corresponding 
ACDM model could be due to systematics. The transition from an accelerating to a decelerating 
expansion occurs at a redshift larger than z = 0.35. Interpreting the dark energy as a minimally 
coupled scalar field we also reconstruct the scalar field and its potential. We constrain Qmo using 
the baryon acoustic oscillation peak in the SDSS luminous red galaxy sample and find that the best 



fit is obtained with f2„ 



0.27, in agreement with the CMB data. 



I. INTRODUCTION 

There is now overwhelming evidence from Supernovae 
data as well as CMB measurements 0,01 and obser- 
vations of large scale structure 0, @ which suggest that 
the Universe is at present undergoing a phase of accel- 
erated expansion. These data also provide evidence that 
nearly two thirds of the total density of the Universe is 
in the form of an effective fluid with a negative pressure. 

Determining the underlying mechanism for this late- 
time acceleration constitutes one of the fundamental 
challenges facing cosmology today. A large number of 
models have been proposed to account for this accelera- 
tion. These fall into two main categories: those that in- 
volve the introduction of an exotic matter source, which 
include the cosmological constant and the quintessence 
scalar fields and those that involve changes to the gravita- 
tional sector of GR, either motivated by String/M theory, 
or through ad hoc modifications of the Hilbert action. 

Among these models some of the most commonly stud- 
ied have been the so called 'dark energy' scenarios, which 
model the underlying accelerating agent as an effective 
perfect fluid with a negative equation of state (EOS). In 
addition to their simplicity, one of the main reasons for 
their popularity has been the fact that an effective model 
of this type can mimic a very wide variety of models, 
ranging from the cosmological constant with EOS equal 
to —1 to models with variable EOS such as quintessence 
models and brane inspired models such as that due to 
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Dvali-Gabadadze-Porrati 0, ^ ■ 

Given the multiplicity of such model candidates, an 
urgent question at present is how to distinguish between 
them observationally. Of particular interest is the nature 
of the EOS of dark energy and importantly whether it is 
variable and different from —1. 

An important approach to this problem has been to 
reconstruct the properties of the dark energy, including 
its equation of state, from the Supernovae 0, and more 
recently the baryon acoustic oscillation (BAO) data 0]. 
A large number of attempts have recently been made at 
such reconstructions. These include schemes that rely 
on specific functional ansatzes for the cosmological pa- 
rameters to be reconstructed, such as for example the 
Hubble parameter ,10^. .llj or the deceleration parameter 
[T3 |. Other schemes employ more general parametrised 
forms such as interpolating fits with right behaviours at 
small and large redshifts ITsI l . There are also the so 
called model-independent schemes which attempt to re- 
cover the cosmological parameters directly from the data 
without assuming their forms by using more general sta- 
tistical tools, such as Markov Chain Monte Carlo tech- 
niques or Gaussian Kernels [l7l |. 

In this paper we put forward a new model-independent 
reconstruction scheme which utilises the expected geo- 
metrical features of the luminosity-distance relation while 
fitting the observational data including the recently re- 
leased Supernova Legacy Survey (SNLS) and baryon 
acoustic oscillation data. Using a phase space represen- 
tation we reconstruct the luminosity-distance curves best 
fitting the data and hence the corresponding cosmologi- 
cal parameters including the Hubble parameter and the 
EOS. 

The plan of the paper is as follows. In section UTI we 
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describe our reconstruction scheme. In Section IIIII we 
give a brief account of the data sets used in our recon- 
structions. Section Hvl cont ains our results: in subsection 
A, we present the results of our reconstructions of the 
cosmological parameters, assuming the cold dark matter 
density parameter to be fimo = 0.27. In subsection B 
we give a discussion of the degeneracy of the luminosity 
distance with respect to flmo ^-nd use the baryon acous- 



tic oscillation data to constrain ri„ 



In subsection C 



we discuss the robustness of our results, and finally we 
conclude in section IVI 



II. THE RECONSTRUCTION METHOD 

A natural starting point in reconstruction schemes 
employing Supernova observations is the luminosity- 
distance di (z) expressed in terms of the Hubble parame- 
ter H: 



di{z) = c{l + z) 



dz 

n ■ 



(1) 



Such reconstructions have taken a variety of forms, often 
relying on specific functional ansatzes for various cos- 
mological parameters. Here we propose an alternative 
method which instead of specifying precise functional 
ansatz for the luminosity-distance relation di{z), imposes 
weak observationally motivated constraints on its geom- 
etry given by: 

(a) d[>0, 

(b) <>0, 

(c) d'r<0, 

where a prime denotes a derivative with respect to the 
redshift z. These are rather natural and weak constraints 
motivated by current observations. In particular condi- 
tion (a) is satisfied for any expanding Universe (dj < 
always implies a contracting Universe). Condition (b) 
is satisfied by any Universe that is currently acceler- 
ating (with a negative deceleration parameter g < 0) 
and which in the past tends to an Einstein-de Sitter 
model with q = 1/2. This can be seen by noting that 
using dH/dt — —77^(1 -I- q), the positivity of c?" im- 
plies q < 1. Finally condition (c) is satisfied by both 
Einstein-de Sitter and KCDA4 models, for all redshifts 
z. This can be seen by employing the statefinder pa- 
rameter r = d/{aH^) 10], and using it to express the 
inequality (c) as r > 3q^ + q — I. This inequality is 
satisfied, for all z, by the AC DM model (for which 
using the Hubble function — poa^^ + A we have 

q=-l + 24poe3^(*+*''V(e^'^^*+*°^ + 4po)^ and r = 1, 
where po and correspond to values at present time). 
The condition saturates for de Sitter space asymptoti- 
cally and is satisfied by the Einstein-de Sitter model (for 



which we have q = 1/2 and r = 1). In terms of the decel- 
eration parameter, condition (c) imposes a lower bound 

^ — 1 

q' on q such that q' > . 

Thus these conditions are compatible with current ob- 
servations and are satisfied, for all z, by both Einstein-de 
Sitter and AC DM models which are the most commonly 
accepted models representing the early and late dynamics 
of the Universe. Together these features provide justi- 
fication for their use in constraining the reconstructed 
luminosity-distance curves, which is the approach we 
shall take in the following. 

Before we proceed with our reconstructions we recall 
that even though the reconstruction process itself does 
not require a theoretical framework, the physical inter- 
pretation of the reconstructed cosmological parameters 
does so. Given that the true underlying theory is not 
known a priori, it is important that such a framework 
is general enough to accommodate a large enough range 
of possible candidates and yet at the same time is suffi- 
ciently simple to make it convenient to work with. Here 
as our theoretical framework we shall adopt general rel- 
ativity plus a minimally coupled scalar field (f> with a 
potential V{(l)). This simple framework includes an im- 
portant set of candidates including GR plus a cosmo- 
logical constant and quintessence models. For simplicity 
we shall consider a flat Friedmann-Lemaitre-Robertson- 
Walker metric. The dark energy can then be viewed as 
a perfect fluid with the corresponding density, pressure 
and EOS given by 



V, P<l> = 



V, w^ = ^, (2) 
P0 



where a dot denotes a derivative with respect to the 
proper time t. 

To facilitate the interpretation of our results, we define 
the following well behaved expansion-normalised vari- 
ables often used in the dynamical systems analysis of 
such models 



(3) 



V 



V6H' 



where p is the matter (cold dark matter plus baryon) 
density. In terms of these variables, the Friedmann and 
the H equations can be written as 



U + V + w — 1 , 

H 3 2 q 

IP ^ r + 

with the equation of state taking the form 

(w^ — v'^) 



W0 



(w^ + v^) 



(4) 
(5) 

(6) 



The reconstruction then proceeds by generating a large 
discrete set, Sdi^ of luminosity-distance curves whose ge- 
ometries are free apart from satisfying the weak con- 
straints (a)-(c) given above as well as fitting the SNLS 



3 



supernovae data set (see the next Section for the details). 
More precisely, each luminosity-distance curve di £ Sdi 
is constructed by a discrete set of points Pi — {zi,di.) 
satisfying the constraints (a)-(c). Let us assume for 
simplicity that Pi are equally spaced points, such that 
dzi — Zi+i — Zi is a constant interval for all i. Then 
each successive point Pi of a luminosity-distance curve 
di € Sdi must satisfy the constraints (a)-(c), which when 
discretised are respectively given by 

di,^^ - 2d;, ^2 -I- di^_^^ < di^^^ - 2di.^^ + di^ (7) 

Thus the reconstruction consists of finding all the 
luminosity-distance curves defined by the discrete set of 
points Pi respecting the above rules and fitting the super- 
novae data. In practice we divided the interval z € [0, 1] 
into steps of dzi = 0.07 for small z and larger steps for 
larger z since the slope of di varies faster for the smaller 
z. 

To reconstruct the discrete set of luminosity distance 
curves, we used an iterative scheme. The starting point is 
a luminosity distance curve defined by the points Pi and 
located in the neighbourhood of and below the SNSL 
data. Then, following the above rules, the iterative 
scheme results in the next luminosity distance curve by 
moving up one of the points Pi. Hence, step by step, the 
slope of the reconstructed curves increase. Each curve 
thus defined is then transformed, using an interpolation 
to go from the discrete form in terms of Pi to a continu- 
ous one. We then check if the reconstructed curve fits the 
data. If so, it is added to the set Sdi- The iterations stop 
when the slope of the reconstructed curve is sufficiently 
large to place it everywhere above the data points, so 
that they are no longer fitted. Note that in the figures 
presented below, interpolation is used to give the recon- 
structed values at equally spaced steps of dzi = 0.07. 
Once the set of curves Sdi is determined, we calculate for 
each curve 



H 



d', 



1 



(8) 



and 



,n I .. 2d, + {l + z)[-2d[ + {l + z)d'^] 

H + {di-{i + z)d[Y • 

Expression (jHJ together with Q can then be used to 
reconstruct u{z) once fimo is specified. Similarly the use 
of expression ® together with the relation d/dt = —(1 + 

z)Hd/dz allows to be reconstructed. Equations Q 
and ^ can then be used to reconstruct the variables 
V and w and hence from ^ the EOS w^. Thus the 
use of this scheme allows the cosmological parameters to 
be reconstructed, without assuming specific ansatzes but 
only the well defined cosmological properties encoded in 
conditions (a)-(c). 



A similar but technically different approach was used 
in 22]. These authors reconstruct the coordinate dis- 
tance which is related to luminosity-distance di by 
the relation di = Hq^hII + z), from supernovae and ra- 
dio galaxy data. In this scheme y is fitted locally, over 
redshift windows of typical length Az k, 0.4, by using a 
quadratic fit to reconstruct the curves y{z) i21j. The first 
and second derivatives of y are then used to reconstruct 
the potential and kinetic energy density. As noted in 
|2ll |. a simpler approach would consist of fitting a poly- 
nomial to the entire set of data. However, low order 
polynomials would not be fiexible enough whereas higher 
order polynomials could induce unphysical oscillatory be- 
haviours. The main difference between our scheme and 
the one used by these authors is that we fit the entire data 
set by imposing a number of geometrical constraints on 
the di curves (and hence implicitly on y) which is not 
necessarily representable by a polynomial in z. This al- 
lows sufficient flexibility in order to fit the entire data 
set while providing enough constraints to eliminate un- 
physical oscillations of di . As we shall see though similar 
our results do not totally agree. For example we obtain 
a much weaker deviation from the ACDM model. 



III. DATA 

In this section we briefly describe the cosmological data 
sets that were used in our reconstructions in the next 
section. 

1. As our primary supernovae sample we took the first 
year data set from the Supernova Legacy Survey 
(SNLS) with 71 new supernovae below z — 1.01, 
together with another 44 low z supernovae already 
available, i.e. a total of 115 SNe. 

2. To check the ability of our scheme to reconstruct 
data, we employed a set of mock data. This consists 
of 115 supernovae with the same distribution in red- 
shift and error bars as the SNLS data but with the 
luminosity-distances replaced by those correspond- 
ing to the ACDM model plus a gaussian noise. 

3. To constrain Qmo we also employed the baryon 
acoustic oscillation peak (BAO) recently detected 
in the correlation function of luminous redgalaxies 
(LRG) in the Sloan Digital Sky Survey 0. This 
peak corresponds to the first acoustic peak at re- 
combination and is determined by the sound hori- 
zon. The observed scale of the peak effectively con- 
strains the quantity 



Ao.35 = D„(0.35)^4r7r^ = 0-469 ± 0.017 , 
0.35c 

where z = 0.35 is the typical LRG redshift and 
Dy(z) is the comoving angular diameter distance 
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defined as 

Dy{z) 

with 



Dm = c 



cz 



dz 



1/3 



H- 



To test the goodness of fit of our reconstructions we em- 
ploy the standard minimization. For the supernovae 
data, is defined by 



E 



where n is the number of data, m°^* and m*'* are respec- 
tively the observed and the theoretically reconstructed 
magnitudes and Oi is the uncertainty in the individual 



IV. RESULTS 

To generate the set of luminosity-distance curves, Sdl^ 
we need to choose the number of points Pi in order to 
describe the curves. We found 7 points to be sufficient for 
the results represented here. See, however, the subsection 
IIV CI below for further discussion of this choice. 

In our reconstructions, each luminosity-distance curve 
was fitted to the SNLS data. Among the reconstructed 
curves Sdi. we only considered those with < 136.8. 
This upper bound on was chosen since for a model 
with one free parameter, such as a flat kCDM model, 
it corresponds to a reasonable value of per degree of 
freedom of Xdof — ^■■^^ .^]- This would clearly be 
larger if > 136.8. 

In the following we shall proceed in two steps. In sub- 
section A we shall assume flmo = 0.27 and reconstruct 
the cosmological parameters using the weak constraints 
(a)-(c) given above. In the subsequent subsection B we 
use the BAO data to constrain il.„ig and show that 0.27 
is indeed the best fitting value. Finally in subsection C 
we study the robustness of our results with respect to the 
number of points used to define the luminosity-distance 
curves, as well as the supernovae data set employed and 
the precise value of flmo used. 



priors on Hq and Hq when making reconstructions using 
only the supernovae data, without assuming additional 
information. 

Here we shall make the following reasonable assump- 
tions. For the present value of the Hubble parame- 
ter we assume 60 < Ho < 80, which from ((SJ im- 
plies that initially 3750 < = ^ < 5000 Mpc. 
Since the derivative of the Hubble parameter is given by 
Hq = ~c{—2d[^ ~^^io ) / ^Zo ^lo ^^^^ small val- 

ues, Hq can take large negative values at present. From 
H'H{1 + z) ^ l/2(/)2 + l/2p, this would correspond to 
high negative values for the EOS Thus in order to 
have a reasonable lower bound on the present value of 
the EOS, we choose Hq > —40, which from ^ implies 
9375 < d['^ = 2dl - H^d'il/c < 13333 Mpc or using 
|3J) corresponds to w^^ > — 2 for the present value of the 
EOS. 

The results presented below were obtained by recon- 
structing more than 23000 luminosity-distance curves fit- 
ting the data. To show that all the reconstructed curves 
respect the Friedmann constraint igj, we have depicted in 
figure^a plot of v'^ +w'^ versus u for these curves, which 
as can be seen provides a very good confirmation of the 
constraint u + v'^ + w'^ = 1. Note that some reconstruc- 
tions with negative energy density {v'^ + 'w'^ < 0) are also 
in agreement with the data, but these correspond to the 
larger values of x^ • This raises the interesting question of 
physicality of such solutions. We shall not consider this 
issue further here, but recall that negative energy densi- 
ties sometimes arise, for instance in quantum cosmology 
or in considerations of braneworld models whose bulks 
possess a negative cosmological constant. 
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A. Reconstruction assuming Qmo = 0.27 

The supernovae data alone do not constrain the present 
values of the Hubble function Hq and its derivative Hq. 
This is due to the fact that there is no data available, and 
hence no constraints, from the future. Consequently, if 
no further information is used, the luminosity-distance 
curves can have any slope at present, subject to the con- 
straints (a)-(c). Thus one needs to assume reasonable 



Figure 1: Plot depicting -t- versus u for the reconstructed 
curves. Tlie perfect straight line provides a good confirmation for the 
constraint u + + — 1. The reconstructed points are shaded 
{coloured: references to colour refer to the web version with colour 
figures) according to their corresponding values, with the deepest 
shade of grey (red), representing the smallest values of ■ 

In figure |2] we have plotted the reconstructed 
luminosity-distance curves (top panel) and the corre- 
sponding Hubble parameters (bottom panel) together 
with their la confidence levels (dashed lines) respectively. 
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As can be seen from the top panel, the di curves best fit- 
ting the data have a slightly smaller (with the small- 
est having = 113.85) than the KCDM model (with 
^ 114). Hence, none of the reconstructed models rules 
out the KCDM model at la. The difference between 
our best fitting luminosity-distance curve and that of 
the ACDM model becomes noticeable at large redshifts 
where the former curve is slightly lower. This deviation 
could come from systematics such as the Malmquist bias. 
Testing our reconstruction method with several samples 
of mock data, we obtain curves which very closely agree 
with the ACDM model even at large redshifts. How- 
ever, they are very similar to the curves obtained using 
the real data. This is due to the fact that the real data 
are themselves very close to the ACDM model, despite 
very small deviations for large redshifts. An example of 
the luminosity-distance reconstruction with mock data is 
displayed in figure O We do not show the plots for the 
other physical quantities, including the EOS, since they 
look very similar to the plots with real data. 

We also find the reconstructed Hubble parameters 
(bottom panel) to have amplitudes which are in quali- 
tative agreement with those found by other authors us- 
ing different reconstruction techniques [T^ . Here again, 
mock data produces similar results. 

Taking as our theoretical framework GR plus a mini- 
mally coupled scalar field, we also reconstructed the ki- 
netic and potential energies corresponding to the scalar 
field (j>, together with their la confidence levels psj . 
These are depicted in figure 0] for the reconstructed di . 
However, to be able to interpret these reconstructions 
within the framework of a minimally coupled scalar field, 
a number of physical constraints need to be taken into 
account jlJl, such as those concerning the signs of the 
potential and kinetic terms. To be cautions, we have ex- 
tracted from the reconstructions the curves correspond- 
ing to the quintessence scalar field models with a positive 
potential. These are depicted in figure[Slshowing that the 
best fitting curves, in the context of quintessence scalar 
field interpretation, are very close to the ACDM model. 
This is particularly striking for z < 0.6 where the degen- 
eracy is very weak. For larger redshifts, some degeneracy 
appears due to larger error bars but remains relatively 
weak. 

We also reconstructed the dark energy EOS (for all the 
reconstructed luminosity-distance curves) whose evolu- 
tion with redshift is depicted in figure To obtain a 
compact representation of the EOS, we found it more 
convenient to have a (u^, w'^) plane representation instead 
of the usual plot of EOS as a function of z. This choice 
of variables is well suited to the particular form 10 of 
the EOS in this case, whose denominator can take small 
values (for instance when the kinetic and potential terms 
take similar values but with different signs) resulting in 
large variations in the EOS as the redshift increases. As 
can be seen from figure |3 this allows a compact repre- 
sentation of the EOS, with the diverging values of EOS 
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Figure 2: Plots of the reconstructed luminosity-distance (top panel) 
and Hubble function (bottom panel) with respect to redshift. The dot- 
ted white (solid black) line represents the ACDM model and the white 
line, the reconstructed luminosity-distance best fitting the data. For 
the Hubble function we have, for sake of clarity, plotted the points 
obtained numerically without joining them. 



corresponding to the dashed line. In this figure the hor- 
izontal line (w-^ = 0) represents the EOS for the cosmo- 
logical constant (w^ = — 1) and the vertical line (w^ = 0) 
the EOS for a stiff fluid (w^ = 1). The remaining line 
defines the w^ = — i/3 line which demarcates the limit 
between the accelerated and decelerated expansion. 

The shaded (coloured) patches in the panels of fig- 
ure IHl represent, from top to bottom, the reconstructed 
values of the EOS at increasing values of the redshift 
given by z = 0,0.28,0.63,0.98. Each point constitut- 
ing the shaded (coloured) patches corresponds to a re- 
constructed luminosity-distance curve, with the deepest 
shades of grey (red) representing the best fits (lowest 
values) . The top panel depicts a cloud of initial points in 
the (u^jW^) plane, which since we are initially assuming 
Mo ~ flmo = 0.27, are forced to lie on a straight line by 
the Friedmann constraint u + v'^ + w'^ — 1. As we go 
to lower panels (higher redshifts), flm evolves differently 
for each point, resulting in the dispersal of the initial 
straight line into different clouds of points. Now since 
we are assuming dj" < 0, c?" becomes smaller as the red- 




Figure 3; Plot depicting tile reconstruction of the luminosity-distance 
relations using the mock data. 



2 . 6 X' 10-"' 
„ 2 . 10-'" 

1 X 10-'" 
^ 5 . 10-1" 



-S xlO-'" 
-1 X 10-1" 





''I 




• 1 . 



0.2 0.4 0.6 0.8 1 



3 ■ 125 . 3,^^^Bll3 . 85 



Figure 4: Plots of all the reconstructed scalar field kinetic (top panel) 
and potential (bottom panel) energies versus the redshift. The horizon- 
tal black lines represent the AC DM model. 



shift increases, but always stays positive. As a result the 
initial straight line configuration of points spreads and 
eventually gets attracted to the neighbourhood of the line 
if^ = + 1/3 which corresponds to d'/ — 0. To see this, 
recah that using (0 and (QJ we have ^d" 2 + HH~^. 
This together with and Q then shows that d'/ = 



implies w 



1/3. As can be seen from the bottom 



panel, this line acts as the accumulation end state of the 
initial cloud of points. 

To summarise, the evolution of EOS with z depicted in 
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Figure 5: Plots of the reconstructed positive scalar field kinetic (top 
panel) and positive potential (bottom panel) energies versus the red- 
shift. The horizontal black lines represent the ACDAl model. 



these panels demonstrates that as z increases, the points 
corresponding to the di best fitting the data stay well 
centred around the cosmological constant line (w^ — 0) 
in agreement with the behaviour found for the scalar field 
kinetic and potential energies. 

We have also depicted in figure the superposition 
of panels similar to those depicted in figure for inter- 
mediate values of the redshift considered in the range 
z G [0, 1], in steps of z = 0.07. We note that some pos- 
sible variation of dark energy has been reported using 
reconstruction techniques different from that employed 
in this paper 0, 0, . Such variations are also pos- 
sible according to our results at la confidence level, but 
are severely constrained in the context of quintessence 
scalar field models, as can be seen from Fig. [S| 

We have also plotted in (top panel of) figure |S1 the 
best fitting reconstructions for the EOS as a function of 
redshift with the la confidence level. As can be seen, ini- 
tially (z ~ 0) the EOS takes values around ~ — 1. Then, 
around the redshift of 0.45, the degeneracy increases be- 
coming very large beyond z = 0.8. This follows the evo- 
lution of error bars and noise levels in the data. We have 
also plotted the evolution of the associated deceleration 
parameter (depicted in the bottom panel of the figure 
|SJ). This shows that a transition from an accelerated to 
a decelerated expansion would occur after the redshift 
of 0.35. There can, nevertheless, exist solutions which 
do not undergo such a transition, even at la confidence 
level, in agreement with [l2l|. This shows that constant 
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Fi gure 6; Figures demonstrating the evolution of the reconstructed 
EOS |6j versus the redshift in the phase plane (u^, w'^). The top panel, 
corresponding to redshift 0, shows the initial points which are forced 
into a linear configuration due to the choice of uq — Ottiq — 0.27 and the 
constraint equation u -\- v'^ -\- nP" — 1. As the redshift increases (lower 
panels), the initial straight line configuration of points spreads and 
eventually gets attracted to the neighbourhood of the line w — -u^ + 1/3 
which corresponds to d[' — 0. 
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Fi gure 7l Figure showing the superposition of similar panels to figure 
l^for all intermediate values of z G [0, 1] in steps of dz — 0.07. As can be 
seen the distribution of the points with the deepest shade of grey (red), 
corresponding to the best fits, is well centred around the cosmological 
constant line. 
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Figure 8: The evolutions of the reconstructed EOS corresponding to 
the luminosity-distance curves best fitting the data as a function of 
the redshift (top panel) and the corresponding deceleration parameter 
(bottom panel) with their la confidence levels (dashed lines). 



redshift z > 0.65. 



as well as negative deceleration parameters are also com- 
patible with the supernovae data, thus indicating that 
at present the supernova data does not establish with 
certainty a transition from accelerated to decelerated ex- 
pansion. 

The panels in figure |51 also show the reconstructed EOS 
and the associated deceleration parameter when the re- 
constructions were confined to the quintessence scalar 
field models with positive potentials. As can be seen at 
la confidence level, the reconstructed EOS is very close 
to that of the KCDM model but could vary weakly be- 
yond z ~ 0.6. Moreover the deceleration can begin for a 



B. Constraining Q.mQ 

Our reconstructions in the previous subsection did not 
provide any information about the CDM density param- 
eter Vtma which was assumed to be 0.27. This is due to 
the fact that the luminosity-distance relation is deter- 
mined purely by the Friedmann equation which is highly 
degenerate with respect to this parameter. To see this, 
let us consider a dark energy model with a constant EOS, 
W0 = r — 1 , and the corresponding Friedmann equation 

{H/H^f = r!™„(l + zf + + zf^ 
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Fi gXire 9: The evolutions of the reconstructed EOS corresponding 
to the luminosity-distance curves best fitting the data as a function of 
the redshift (top panel) and the corresponding deceleration parameter 
{bottom panel) with their la confidence levels (dashed lines), using the 
scalar field representation confined to quintessence models. 



Re- writing the CDM density parameter as Qr 
^imo + ^2mo5 the Friedmann equation becomes 



ar 



(H/Hof = Quno (1 + + ^2mo (1 + + Q^, (1 + z) 

This form of H{z) may be viewed as a new dark energy 
model with a different CDM density parameter Slim^, 
and a different dark energy density represented by 
the last two terms in this expression. Hence without 
changing the Hubble function we can re-group the last 
two terms as a new dark energy term thus: 

P<Po 

with a corresponding equation of state given by 

^^^^^ n2m„{i + zf + n^,{i + zr^ 

Thus given a Hubble function we can construct different 
representations with different effective CDM density pa- 
rameters and dark energy components, but with identical 
luminosity-distance functions di. This makes transpar- 
ent the fact that the luminosity-distance di is potentially 
highly degenerate with respect to the CDM density in 
Universe. It also illustrates the fact that a constant EOS 
can misleadingly become time-dependent if the matter 
density is incorrectly chosen, as is shown in p^ . We 
should, however, add that since in practice ftmo cannot 
be determined precisely, this demonstrates that a real 
observer can never prove that w^ is truly constant, simi- 
lar to the limitations that exist in determining the exact 
flatness of the Universe from observations. In practice 
one can only take the most likely value for flmo given 



by the most accurate available observations in order to 
determine the corresponding EOS. 

To constrain r^mo, we need further input from ob- 
servations. Here we used the recent baryon oscillation 
data loj (discussed in Section III), to further constrain 
each luminosity-distance curve. We found that for our 
reconstructed set of curves Sdh lies in the range 

0.16 < ilmo < 0.41, with the best fit obtained for 
^nio = 0.27. This can be seen clearly from the plot 
of flmg versus shown in figure ^1 which justifies the 
choice of the specific value ftmo = 0.27 in the previous 
subsection. 
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Figure 10: Figure depicting the variation of reconstructed with 
respect to i using tlie SNLS and BAO data. 



C. Robustness of the reconstruction 

In this subsection we study the robustness of our re- 
sults obtained in the previous subsections. 

In our reconstructions so far we have used 7 points 
{di,z) to describe the luminosity-distance curves. This 
choice was made by trial and error. We find that by tak- 
ing too many points we can "over-fit" the data, in the 
sense that we can always obtain a perfect but artificial 
match to any set of data. In that case one would in effect 
be fitting noise with an anomalously low best which 
can, among other things, rule out the KCDM model at 
Icr. On the other hand, taking too few points would 
amount to " under- fitting" the data. Thus one would not 
be able to fit the model represented by the data and the 
best can be anomalously high. For example, it could 
make it impossible to fit the KCDM model. In order 
to determine the optimal number of points, the use of 
mock data is extremely useful. We have tested our re- 
construction scheme with several sets of mock data, all 
based on the KCDM model. Using seven points we were 
able to recover, in all cases, the KCDA4 model at la 
with a best (typically around 113.7) only a few tenths 
of percent smaller than 114 (which is close to the cor- 
responding value for the KCDM model). This demon- 
strates that, given the distribution of the SNLS data and 
the corresponding error bars, our reconstruction process 
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works well (although a slight over-fitting is unavoidable 
in practice as one is not expected to recover the best 
fit value exactly). We checked that our results remain 
robust with small changes in the number of points. 

To summarise, care must be taken in choosing the 
number of points taken to define a luminosity-distance 
curve. We found that the employment of the mock data 
together with our conditions (a)-(c) allow an appropriate 
number of points to be chosen, i.e. around 7. 

We note that the number of points defining a 
luminosity-distance curve is not related to the degrees 
of freedom of the underlying theory giving rise to that 
curve. Thus, a straight line can be defined by 7 points 
even though only 2 points are necessary to construct the 
line, while analytically the equation of a straight line 
passing through the origin only has one degree of free- 
dom. 

We also checked that reasonable changes in the value 
of does not alter the qualitative behaviour of dark 
energy EOS found here. 

V. CONCLUSION 

We have proposed a model-independent reconstruc- 
tion scheme which is compatible with current observa- 
tions and shares a number of geometrical features with 
the Einstein-de Sitter and ACDM models, which are the 
most commonly accepted models representing the early 
and late dynamics of the Universe. Together these fea- 
tures provide justification for our proposed scheme. 

Using this scheme, and assuming ilmo = 0.27 together 
with the SNLS supernovae data, we have reconstructed 
a large set of luminosity-distance curves. Using these 
curves we have reconstructed the cosmological parame- 
ters, including the EOS. Our reconstructions show that 
the luminosity-distance curves best fitting the data cor- 
respond to a slightly varying dark energy density with 
the Universe expanding slightly slower than the ACDM 
model. However, the ACDM model fits the data at 
1(7 and the fact that our best fitting luminosity-distance 
curve is lower than that of the corresponding ACDM 
model could be due to systematics such as Malmquist 
bias. Reconstructing the EOS, large degeneracy appears 
around z — 0.6 in agreement with increasing error bars 
and noise levels in the distribution of data. The recon- 



structed deceleration parameter shows that the transition 
redshift to a decelerating Universe should be larger than 
z = 0.35. 

Assuming the theoretical framework to be GR plus a 
minimally coupled scalar field, we also considered recon- 
structions confined to quintessence models with positive 
potentials. In that case we find the best fitting recon- 
structed curves to be very close to those corresponding to 
the ACDM model, in particular for low redshifts z < 0.6. 
For larger redshifts, larger error bars in the data allow a 
small increase in the kinetic term and EOS and a small 
decrease in the potential term at la. However, in the con- 
text of quintessence interpretation, it seems rather hard 
to find a model fitting the data better than the ACDM 
model at la. 

We also used the BAO data to constrain ilmo and 
found that takes values in the range 0.16 < ftm„ < 
0.41 with a best fit given by il.mg = 0.27, in close agree- 
ment with CMB data 

Using different techniques, and emplo ying supernovae 
as well as other data, other authors 0i ll^l '2^ also 
find evidence for a possible variation of dark energy den- 
sity. However, a constant dark energy density is again 
never ruled out to a satisfactory confidence level and thus 
the important question of whether the EOS varies in time 
or not still remains open. 

The present paper takes a new step in reconstruct- 
ing the equation of state by providing a new model- 
independent reconstruction scheme of the cosmological 
parameters. Our approach can be applied to other the- 
oretical frameworks which can be cast into GR plus a 
perfect fiuid. It would be interesting to repeat these 
reconstructions with the next set of SNLS data to de- 
termine whether the best fitting curves move closer or 
further away from the ACDM model. 
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